Quasars with a Kick — Black Hole Recoil in Quasars 



G. A. Shields 2 , E. W. Bonning 3 , and S. Salviander 2 



ABSTRACT 

Mergers of spinning black holes can give recoil velocities from gravitational radiation 
up to several thousand km s _1 . A recoiling supermassive black hole in an AGN can 
retain the inner part of its accretion disk, providing fuel for continuing AGN activity. 
Using AGN in the Sloan Digital Sky Survey (SDSS) that show velocity shifts of the 
broad emission lines relative to the narrow lines, we place upper limits on the incidence 
of high velocity recoils in AGN. Brief but powerful flares in soft X-rays may occur when 
bound material falls back into the moving accretion disk. 

Subject headings: galaxies: active — quasars: general — black hole physics 



1. Introduction 

Recent numerical simulations of binary 
black hole mergers by a number of groups 
show large recoil velocities or "kicks" result- 
ing from from anisotropic emission of gravi- 



tation al radiation [see, e.g.. ICampanelli et al 



(120071 ) for a summary and references]. These 
results show a maximum kick for spins anti- 
aligned and perpendicular to the orbital an- 
gular momentum, re aching 2500 km s -1 
for spin a* = 0.8 ([Gonzalez et al.1 12 007 ; 



Tichy fc Marronettill2007l ). ICampanelli et al 
(120071 ) predict a maximum recoil velocity of 
4000 km s _1 for equal mass black holes with 
maximal spin. For a binary supermassive 
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black hole ( ~ 10 8 Mm) formed duri ng a galac- 



tic merger (|Begelman et al.1 ll98Ch . the kick 



may displace the black hole fr om the galac- 



tic n ucleus or eject it entirely (jMerritt et al 
20041 . and references therein). 



If an accretion disk is present, the inner 
part will remain bound to the black hole, fu- 
eling AGN activity as the black hole leaves the 
galactic nucleus . For a merger taking place in 
an AGN, the accretion disk fueling the AGN 
will remain bound to the recoiling black hole 
inside the radius = 10 18,1 Mg/f 2 000 cm 
where the orbital velocity equals fkkk • Here, 
M s = M BH /1O 8 M and v mm = vl^/WOO 
km s -1 . For an q-disk (jShakura fc Sunvaev 



1973 : iFrank et al J 12002 ). the disk mass M b 
that remains bound is 



M b 



(10 8 -°M o )a:y 5 M a 3/2 M n 7/10 ^ 5/2 



T000 



(1) 



where Mo is the accret i on rate in solar masse s 
per year (jLoebl 120071 : iBonning et all 120071 ). 
subject to M(j is k < Mbh for stability. Until 
the disk is depleted, the post-merger accre- 
tion rate will resemble the pre-merger accre- 
tion rate, giving a disk consumption time td ~ 



1 



M b /M « (10 8 yrJal^Aff^Mo- 3 / 10 ^ 2 . 
Thus a recoiling AGN could shine for a time 
comparable with estimates of the total AGN 
lifetime. Here we discuss two observational 
manifestations of black hole recoils in AGN, 
involving Doppler shifts of the emission lines 
and energy released by infall into the disk fol- 
lowing the recoil. 

2. Shifted emission lines 

The broad emission-line region (BLR) of 
AGN corresponds to material that remains 
bound to the black hole, whereas the nar- 
row lines arise from gas in the potential of 
the host galaxy that will not follow the recoil- 

" 20061 ). There- 



ing black hole (jMerritt et al 



fore, the broad lines will appear shifted with 
respect to the galaxy systemic redshift as ex- 
pressed by the narrow emission lines. We 
have carried out a search for candidate kicked 
QSOs using spectra from the Sloan Digital 
Sky Surve>a Data Release 5 (DR5). We used 
objects spectroscopically identified as QSOs 
by SDSS having 0.1 < z < 0.81 such that H/3 
and [O in] A5007 were both measurable. The 



lines were measured as in ISalviander et al 
(120071 ) . involving subtraction of Fe II and nar- 
row H/3, and a Gauss-Hermite fit to the line 
profiles with asymmetry parameter /13 and 
kurtosis parameter /14. Our final data set con- 
sists of 2598 obj ects. This reflect s the v arious 
quality cuts of Salviander et al. ( 20071 ) plus 
elimination of objects with strongly asymmet- 
rical H/3 giving /13 greater than 0.1 in absolute 
value. The displacement of the broad H/3 line 
peak with respect to the peak of [O in] was 
calculated as Avjip = c(znp — zq m)/(l + 
zp ttt) and analogously for Mg II. For details 
see 



Bonning et all (|2007h . 



The distribution of Auwg has a mean dis- 
placement of +100 km s _1 and a FWHM of 
500 km s _1 . The distribution of AvMg has 



a similar FWHM and mean displacement of 
-150 km s _1 . In our data set of 2598 objects, 
the fractions f v having Auh/3 greater than a 
given value are: / 50 o = 0.04, / 800 = 0.007, 
/1000 = 0.0035, /1500 = 0.0012, / 200 o = 
0.0008, and / 250 o = 0.0004. However, we re- 
gard these values as upper limits on the num- 
ber of true recoils: (1) Mg 11 typically has 
a smaller shift than H/3, such that A^Mg ~ 
0.6Avh/3 • For recoils, one expects similar 
shifts for all broad lines. (2) The largest shifts 
occur only for objects with large H/3 FWHM, 
a correlation not expected for recoils. (3) The 
high Atjh/3 objects generally have normal nar- 
row line intensities and [Ne v] profiles wider 
than the lower ionization lines, consistent with 
a normal NLR with a central ionizing source. 
(4) Some high Avup objects show a charac- 
teristic asymmetry involving a suppression of 
the flux on one side of the line, suggesting 
complications in the BLR. 

Among the 501 objects having measured 
velocities for both H/3 and Mg II, there are no 
objects with AvMg > 800 km s _1 for which 
Avh/3 and Ai>M g agree within expectation. 
This suggests /soo < 0.002, a tighter limit 
than above. 

For mergers of two black holes with spin 

a* = 0.9 an d a range of mass ratio q = 

rri2/mi > 0.2, Schnittman Buonanno ( 20071 ) 



estimate theoretical values /, 

/; 



500 



0.31 and 



1000 



0.079. Convolving with random kick 
inclinations, we find the predicted line-of-sight 
fractions to be /500 = 0.18 and /1000 = 0.054. 
Our observational results give upper limits on 
f v substantially lower than these predictions. 
However, a number of factors may contribute 
to the low observed incidence of recoil shifts. 
(1) A thorough exploration of parameter 
space for mergers is needed using full numer- 



ical relativity simulations (jCampanelli et al 



20071 ). (2) The distribution of spin parame- 



J The SDSS website is |http : / /www .sdss.org 



ters a* for supermassive black holes remains 
uncertain, as does the expected distribution of 



2 



merger parameters ( mass ratio, spin orienta - 
tions). For example. iBogdanovic et al.l (|2007l ) 
argue that gas accretion will align the spins 
with the orbital angular momentum vector, 
resulting in small kicks. (3) The likelihood of 
black hole mergers early in the luminous phase 
of an AGN is uncertain. (4) The observed 
shift velocity may be substantially reduced 
by dynamical effects including the inertia of 
the bound disk and the stars that remain 
bound to the recoiling hole, as well as the 
slowing of the black hole by the gravitational 
potential of the galaxy and by dynamical fric- 
tion o yer the luminous life of the wandering 
AGN (jMerritt et al.ll2004l : iMadau & Quataert 
20041 ). 

These dynamical effects are less severe for 
large kick velocities. For our sample, the me- 
dian va lue of Mrh, derived fro m the H/3 line 



width (jSalviander et al.1 120071 ') is 10 Y y M , 
with L/L Ed « 10-°- 8 , and M = lO" 5 . 
Here we have assu med L = O.lMc 2 and 
L bo i = 9i/L„(5100) (jKaspi et alJ l200d ). For 



this configuration, the bound disk mass is 



Mb (10 7 - 4 M Q ) v^qJq. The momentum of 
the recoiling black hole will be shared by the 
mass of the bound disk traveling with it, giv- 
ing v = UkickMBH/(M B H + M b ). However, 
for Wkick > 1500 km s _1 this does not seri- 
ously slow the hole. The typical Mbh for our 
sample corresponds to <r* ~ 175 km s -1 by 
the M bh — c* relationship of iTremaine et al 



h00i ). A recoil at 



^kick 



2000 km s" 



gives Mb = 0.07Mbhj which slows the hole to 
1900 km s^ 1 . The disk consumption time is 
about 17 million years. We have estimated the 
slowing of the hole in a simple model for the 
host gal axy resembling the i s other mal sphere 
model of Madau fc Quataertl (|2004h . The hole 
reaches a radius of roughly 30 kpc while slow- 
ing to ~ 1500 km s _1 during the luminous 
phase. (For this high velocity, dynamical fric- 
tion is minor during td-) Most of the luminous 
phase is spent close to this velocity, and half 



of the objects will have a line-of-sight veloc- 
ity component above 800 km s _1 for random 
orientations. 

If the AGN in our sample have a typical 
lifetime of 50 million years, whereas the life 
of objects with i>ki c k > 2000 km s _1 is about 
17 million years, then the fraction of AGN 
with an observed shift over 800 km s^ 1 will 
be /goo ~ (0.5)(17/50)/| 000 , where /| 000 is 
the intrinsic fraction of recoils over 2000 km 
s _1 . This may be compared with our ob- 
served limit of less than 0.2% for shifts over 
800 km s -1 , implying that f^ 000 < 1.2%. This 
example suggests that the statistics of velocity 
shifts have potential to give useful constraints 
on high velocity recoil events. 

3. Recoil flares 

Marginally bound material will lag behind 
the recoiling hole and then fall back into the 
disk with velocity ~ Wkick- In the frame of 
the moving hole, the bound disk material has 
excess energy £ flare = (l/2)M b t^ ick compared 
to material in circular orbit with the pre-kick 
angular momentum. The energy associated 
with the mass given by Eq. [I] is 

£ flare = (10 57 - 3 erg)alf M^M^v^. 

If this energy is dissipated on the dynamical 
timescale 



tfla 



Rb/v ki ck = (420 yr)M 8Vl - 3 00 (3) 



by shocks as the material falls back, the power 
associated with this "recoil flare" is 



L 



flare 



(10 47 - 2 erg s->:7 u M 8 1/Z M ( 



4/5 ,,1/2 .-,7/10 5/2 







y iooo- 
(4) 



The post-shock temperature will be 

T shock = (10 6 - 9 K)( Ulo0 o) 2 = (0.7 keV) V 2 000 

(5) 

(Osterbrock & Ferland 2006). The flare may 
substantially exceed the luminosity of the 



3 



QSO in its active phase as well as the Edding- 
ton limit, raising the possibility of radiation 
driven outflows. 

The shocks will emit soft X-rays, but ab- 
sorption by the infalling material will be 
severe. Consider a QSO with Mbh = 
10 8 M accreting at lO -05 M yr" 1 ( giving 
L/L^d ~ 10 -0 ' 9 ) before the final black hole 
in-spiral. For a kick velocity of 1000 km s _1 , 
the bound disk has mass lO 7 ' 7 M , radius 
i?b = 10 18 ' 1 cm, and surface density S ~ 
10 4 ' 2 g cm -2 . This corresponds to an electron 
scattering optical depth (if ionized) r e ~ 10 3 ' 8 . 
The optical depth due to photoelectric ab- 
sorption by H, He, and heavy elements will be 
large. The density of the infalling material is 
N = 10 9 ' 9 atoms per cubic centimeter. For a 
rough estimate of the spectrum of the cooling 
shocked gas, we computed a coronal equilib- 
rium model using version 07.02.00 of the pho- 
toionization code CLOUDY, mo st recently 
described by Ferland et al. ( 19981 ). We used 



T = 10 7 K, N = lO 10 6 cm" 3 , and solar abun- 
dances. The cooling time is ~ 10 -3 ' 1 yr, with 
much of the total cooling in the form of ther- 
mal bremsstrahlung. There are many strong 
emission lines, including lines of Fe XVII to 
Fe XX at 12 to 15 A, each with several per- 
cent of the total energy. The thickness of 



10 



11.8 



cm, small 



the cooling zone will be 
compared with R\>. 

The X-ray luminosity will photoionize the 
infalling material approaching the shock. The 
local energy flux in the X-ray flare is Ffiare ~ 
10 9 ' 8 erg s™ 1 cm -2 , giving an ionizing pho- 
ton flux of 4>i ~ 10 19,4 cm -2 . We have com- 
puted photoionization models with CLOUDY 
for this precursor zone. We used a thermal 
bremsstrahlung ionizing spectrum at 10 7 K 
with flux 10 10 erg s~ 4 cm~ 2 , a gas density 
N = 10 10 cnr 3 , and a cutoff column density 
10 24 cm -2 . The photoionized depth in this 
pre-shock flow is ~ 10 13 ' 5 cm, again narrow 
compared to Rf,- The high ionization param- 



eter (U = 10~°' 9 ) and hard ionizing spectrum 
give strong ultraviolet emission lines of highly 
ionized species. The Lya line has an emer- 
gent flux of 8% of the incident flux, 28 times 
the flux of H/3. Other line intensities include 
I(X)/I(Lya) = 1.5 for O VI A1035 and 0.8 
for C IV A 1549. The high density suppresses 
forbidden line emission. The line widths will 
be ~ fkick- Any surrounding dusty torus (An- 
tonucci & Miller 1985) will in turn reprocess 
much of this energy into infrared radiation 
from dust. 

The number of recoil flares currently ob- 
servable involves the rate of mergers, the 
probability of kicks of various velocities, and 
the flare duration at a given velocity. Let 
dN/dt be the merger rate per year of observer 
time in some range of Mbh and redshift z. 
Roughly, the number of kicks currently ob- 
servable is iV flarc Ri (dN/dt) [t flarc (l + Z)]fl, 

where is as above. Haehnelt (1994) derived 
the black hole merger rate from the QSO lu- 
minosity function. For Mbh between 10 7 and 
10 s Mq, Haehnelt finds 0.016 events per year 
of observer time for z < 3, giving roughly 6, 
0.2 flares in play for t^ick > 500, 1000 km s , 
respectively. Uncertainties include the flare 
duration, the number of obscured QSOs, and 
Haehnelt's assumed L/Lfid an d QSO lifetime. 

For a 10 8 M hole in a QSO shining at 
0.3LEd at z = 2, the received flux from the 
flare is ~ 10 -1 erg s -1 cm" 2 for viooo = 1- 
The observed spectrum will be softened by the 
redshift, but a considerable fraction of the ra- 
diation should still be received at photon en- 
ergies above 0.2 keV. At 10 -11 ' 7 erg s _1 cm -2 , 
Hasinger, Miyaji, & Schmidt (2005) find ~ 0.2 
soft X-ray sources per square degree, so that 
the recoil flares would need to be identified 
from among roughly 10 4 other sources of sim- 
ilar flux. 

Recoil flares may provide a unique oppor- 
tunity to detect high velocity recoils in AGN. 
Detailed simulations of the gas dynamics and 



4 



emitted spectrum are needed to help identify 
these rare events. 
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